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Abstract 

The renormalization group method of Goldenfeld, Oono and their collaborators 
is applied to asymptotic analysis of vector fields. The method is formulated on the 
basis of the theory of envelopes, as was done for scalar fields. This formulation 
actually completes the discussion of the previous work for scalar equations. It is 
shown in a generic way that the method applied to equations with a bifurcation 
leads to the Landau-Stuart and the (time-dependent) Ginzburg-Landau equations. 
It is confirmed that this method is actually a powerful theory for the reduction of 
the dynamics as the reductive perturbation method is. Some examples for ordinar 
diferential equations, such as the forced Duffing, the Lotka-Volterra and the Lorenz 
equations, are worked out in this method: The time evolution of the solution of the 
Lotka-Volterra equation is explicitly given, while the center manifolds of the Lorenz 
equation are constructed in a simple way in the RG method. 



1 Introduction 



It is well known that the renormalization group (RG) equations |1[] have a peculiar power 
to improve the global nature of functions obtained in the perturbation theory in quantum 
field theory (QFT)||: The RG equations may be interpreted as representing the fact 
that the physical quantities 0(p,a,fi) should not depend on the renormalization point /i 
having any arbitrary value, 

dOjp, ol\h) = Q 
dfj, 
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Such a floating renormalization point was first introduced by Gell-Mann and Low in the 
celebrated paper Q. 

It is Goldenfeld, Oono and their collaborators ( to be abbreviated to GO) || f| who 
first showed that the RG equation can be used for purely mathematical problems as 
to improving the global nature of the solutions of differential equations obtained in the 
perturbation theory. One might say, however, that their presentation of the method is 
rather heuristic, heavily relied on the RG prescription in QFT and statistical physics; it 
seems that they were not so eager to give a mathematical reasoning to the method so that 
their method may be understandable even for those who are not familiar with the RG.Q 
In fact, the reason why the RG equations even in QFT "improve" naive perturbation had 
not been elucidated. One may say that when GO successfully applied the RG method 
to purely mathematical problems such as solving differential equations, it had shaped a 
clear problem to reveal the mathematical reason of the powefullness of the RG method, 
at least, a la Stuckelberg-Peterman and Gell-Mann-Low. 

Quite recently, the present author has formulated the method and given the reasoning 
of GO's method on the basis of the classical theory of envelopes ||, ||: It was demonstrated 
that owing to the very RG equation, the functions consturcted from the solutions in the 
perturbation theory certainly satisfies the differential equation in question uniformly up to 
the order with which local solutions around t — to is constructed. It was also indicated in 
a generic way that the RG equation may be regarded as the envelope equation. In fact, if 
a family of curves {C^}^ in the x-y plane is represented by y — f(x; fi), the function g(x) 
representing the envelope E is given by eliminating the parameter \i from the equation 



One can readily recognize the similarity of the envelope equation Eq.(1.2) with the RG 
equation Eq.(l.l). In Ref.'s|5], [|, a simplified prescription of the RG method is also 
presented. For instance, the perturbative expansion is made with respect to a small 
parameter and independent functions^, and the procedure of the "renormalization" has 
been shown unnecessary. 

However, the work given in []5], || may be said to be incomplete in the following sense: 
To give the proof mentioned above, the scalar equation in question was converted to a 
system of first order equations, which describe a vetor field. But the theory of envelopes 
for vetor fields, i.e., envelopes of trajectories, has not been presented in |5|, ||. The theory 
should have been formulated for vector equations to make the discussion self-contained 
and complete. 

One of the purposes of the present paper is therefore to reformulate geometrically the 
RG method for vector equations, i.e., systems of ODE's and PDE's and to complete the 

1 In Appendix A, we give a brief account of the Goldenfeld et al's prescription. 

2 Such an asmptotic series is called generalizes asymptotic series. The author is gratcfull to T. Hatusda 
for telling him this fact and making him recognize its significance. 



df(x;fj) 
djj, 



0. 
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discussion given in @, 0. 



Another drawback of the previous work is that a reasoning given to a procedure to 
setting to = t in the RG method^ of Goldenfeld et al was not fully persuasive. f\ In this 
paper, we present a more convincing reasoning for the procedure. 

Once the RG method is formulated for vector fields, the applicability of the RG method 
developed by Goldenfeld, Oono and their collaborators is found to be wider than one might 
have imagined: The RG method is applicable also to, say,n-dimensional vector equations 
that are not simply converted to a scalar equation of the n-th order; needless to say, it 
is not necessarily possible to convert a system of ordinary differential equations (or dy- 
namical system) to a scalar equation of a high order with a simple structure, although 
the converse is always possible trivially. For partial differential equations, it is not always 
possible to convert a system to a scalar equation of a high order 0]. Moreover, interest- 
ing equations in science including physics and applied mathematics are often given as a 
system. Therefore, it is of interest and importance to show that the RG method can be 
extended and applied to vector equations. To demonstrate the powefulness of the method, 
we shall work out some specific examples of vector equations. 

We shall emphasize that the RG method provides a general method for the reduction 
of the dynamics as the reductive perturbation method (abbreviated to RP method) || 
does. It should be mentioned that Chen, Goldenfeld and Oono[|J already indicated that 
it is a rule that the RG equation gives equations for slow motions which the RP method 
may also describe. In this paper, we shall confirm their observation in a most general 
setting for vector equations. Furthermore, one can show [[| that the natural extension of 
the RG method also applies to difference equations or maps, and an extended envelope 
equation leads to a reduction of the dynamics even for discrete maps. Thus one sees that 
the RG method is truly a most promising candidate for a general theory of the reduction 
of the dynamics, although actual computation is often tediuous in such a general and 
mechanical method. 



This paper is organized as follows: In the next section, we desribe the theory of en- 
velopes for curves (or trajectories) in parameter representation. In section 3, the way 
how to construct envelope surfaces is given when a family of surfaces in three-dimensional 
space are parametrized with two parameters. In section 4, we give the basic mathe- 
matical theorem for the RG method applied to vector fields. This section is partially a 
recapitulation of a part of Ref. ||, although some clarifications are made here. In sec- 



tion 5, some examples are examined in this method, such as the forced Dufnng(T(|, the 
Lotka-Volterra||TTJ and the Lorenz [|T^, |I~0|j equations. The Duffing equation is also an ex- 
ample of non-autonomous one, containing an external force. In section 6, we treat generic 
equations with a bifurcation; the Landau-Stuart [TBj and the Ginzburg-Landau equation 
will be derived in the RG method. The final section is devoted to a brief summary and 
concluding remarks. In Appendix A, a critical review of the Goldenfeld et al's method 



3 See Appendix A. 

4 The author is grateful to Y. Oono for his criticism on this point. 
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is given. In Appendix B, the Duffing equation is solved as a scalar equation in the RG 
method. 



2 Envelopes of trajectories 



To give a geometrical meaning to the RG equation for systems, one needs to formulate a 
theory of envelopes of curves which are given in a parameter representation: For example, 
if the equation is for u(t) = t (x(t),y(t)) ) the solution forms a trajectory or curve in 
the x-y plane with t being a parameter. In this section, we give a brief account of the 
classical theory of envelopes for curves in the n-dimensional space, given in a parameter 
representation. 

Let a family of curves {C a } a in an n-dimensional space be given by 

X(t;a)= t (X 1 (t;a),X 2 (t;a),...,X n (t;a)), (2.1) 

where the point (X 1 , X 2 , ■■■X n ) moves in the n-dimensional space when t is varied. Curves 
in the family is parametrized by a. We suppose that the family of curves {C a } a has the 
envelope E: 

X E (t) = t (X E1 (t; a),X E2 (t; a), X En (t; a)). (2.2) 



The functions X E (t) may be obtained from X(t; a) as follows. If the contact point of 
C a and E is given by t = t a , we have 

X(t a ,a) = X E (t a ). (2.3) 

For each point in E, there exists a parameter a = a(t): Thus the envelope function is 
given by 

X E (t a )=X(t a ;a(t a )). (2.4) 



Then the problem is to get the function a(t), which is achieved as follows. The 
condition that E and C a has the common tangent line at X(t a ; a) = X E (t a ) reads 





dX 
~dJ 


t — t a 


dX E 
dt 






On the other hand, differentiating Eq.(2.4), one 


has 




dX E 
dt 


t — tot 


OX 
~dt 


dX da 

t=t a ' da dt 





(2.5) 



(2.6) 
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From the last two equations, we get 

From this equation, the function a = a(t) is obtained. This is of the same form as 
the RG equation. Thus one may call the envelope equation the RG/E equation, too. 
In the application of the envelope theory for constructing global solutions of differential 
equations, the parameter is the initial time t , i.e., a = t . Actually, apart from t , we 
have unknown functions given as initial values in the applications. We use the above 
condition to determine the to dependence of the initial values by imposing that t = t. 
In section 4, we shall show that the resultant function obtained as the envelope of the 
local solutions in the perturbation theory becomes an approximate but uniformly valid 
solution. 



3 Envelope Surfaces 

This section is devoted to give the condition for constructing the envelope surface of a 

family of surfaces with two parameters in the three-dimensional space. The generalization 
to the n-dimensional case is straightforward. 

Let {S rir2 } riT2 be a family of surfaces given by 

F(r;r 1 ,r 2 )=0, (3.1) 
and E the envelope surface of it given by 

G(r) = 0, (3.2) 

with r = (x, y, z). 

The fact that E contacts with S T1T2 at (x, y, z) implies 

G(r) = F(r;r 1 (r),r 2 (r)) = 0. (3.3) 

Let (r + dr, n + dn, t 2 + dr 2 ) gives another point in E, then 

G(r + dr) = F(r + dr; n + dn, r 2 + dr 2 ) = 0. (3.4) 

Taking the difference of the two equations, we have 

dF dF 

VF-dr + —dn + —dr 2 = 0. (3.5) 

On the other hand, the fact that E and S rir2 have a common tangent plane at r implies 
that 

VF-dr = 0. (3.6) 
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Combining the last two equations, one has 



dF dF 

£-dn + ^-dr 2 = 0. (3.7) 
on dr 2 

Since dr\ and dr 2 may be varied independently, we have 

dF dF . . 

— = 0, — = 0. 3.8 
oti dr 2 

From these equations, we get as a function of r; n = Tj(r). 
As an example, let 

F{x, y, z; n, r 2 ) = e~™{l - y(x - n)} + e"^{l - x(y - r 2 )} - z. (3.9) 
The conditions dF ' jdr\ = and dF jdr 2 = give 

n = x, r 2 = y, (3.10) 
respectively. Hence one finds that the envelope is given by 

G{x, y, z) = F(x, y, z; n = x, r 2 = y) = 2e"^ - z = 0, (3.11) 
or z = 2exp(— xy). 

It is obvious that the discussion can be extended to higher dimensional cases. In 
Ref. 0, envelope surfaces were constructed in multi steps when the RG method was applied 
to PDE's. However, as has been shown in this section, the construction can be performed 
by single step. 

4 The basis of the RG method for systems 
4.1 ODE's 

LetX= t (X 1 ,X 2 ,---,X n ) andF(X,t;e) = '(^(X, t; e), F 2 (X,t; e), • • • , F n (X, t; e)), and 
X satisfy the equation 

f = F(X, i;£ ). (4.1) 

Let us try to have the perturbation solution of Eq.(4.1) around t = t by expanding 

X(t; t ) = X„(t; t ) + eXi(t; t ) + e 2 X 2 (t; t ) • • • . (4.2) 
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We suppose that an approximate solution X = X(£; to, W(t )) to the equation up to 0(e p ) 
is obtained, 

^'f°»L F (x W , 1;t ) + 0W . (4.3) 

where the n-dimensional vector W(to) denotes the initial values assigned at the initial 
time t = to. Here notice that to is arbitrary. 

Let us construct the envelope function X^(t) of the family of trajectories given by 
the functions X(t; t , W(t )) with t parameterizing the trajectories. The construction is 
performed as follows: First we impose the RG/E equation, which now reads 

— = 0. (4.4) 
ato 

Notice that X contains the unknown function W(to) of to-Q in the usual theory of en- 
velopes, as given in section 2, this equation gives to as a function of t. However, since we 
are now constructing the perturbation solution that is as close as possible to the exact 
one around t = t , we demand that the RG/E equation should give the solution t = t, 
i.e., the parameter should coincide with the point of tangency. It means that the RG/E 
equation should determine the n-components of the initial vector W(t ) so that t = t. 
In fact, Eq.(4.4) may give equations as many as n which are independent of each other.Q 
Thus the envelope function is given by 

X E (t)=X(t;t,W(t)). (4.5) 



Then the fundamental theorem for the RG method is the following: 
Theorem: X E {t) satisfies the original equation uniformly up to 0(e p ). 

Proof The proof is already given in Eq.(3.21) of Ref. ||. Here we recapitulate it for 
completeness. Vto, owing to the RG/E equation one has 



dX f 



dt 



dX(t;t ,W(t )) 



t=t 



dt 

dX(t;to,W(t )) 



dX(t;t ,W(t )) 



t=t 



dt 



t=t 



dt r j 

= F(X E (to),t ;e) + 0(e p ), 

where Eq.(4.4) has been used in the last equality. This concludes the proof. 
5 This means that Eq.(4.4) is a total derivative w.r.t. to; 

d± d± dw ax 



dt dt dt dW 



(4.6) 



3 In the applications given below, the equation is, however, reduced to a scalar equation. 
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4.2 PDE's 



It is desirable to develop a general theory for systems of PDE's as has been done for 
ODE's. But such a general theorem is not available yet. Nevertheless it is known that 
the simple generalization of Eq. (4.4) to envelope surfaces works. 

Let X(£, x : to, x o; W(to> x o)) is an approximate solution given in the perturbation 
theory up to 0(e p ) of a system of PDE's with respect to t and x = (x%, x%, . . . , x n ). Here 
we have made explicit that the solution has an initial and boundary value W(to,xo) 
dependent on t and x = (x 10 , x 20 , . . . , x n0 ). As has been shown in section 3, the RG/E 
equation now reads 

^ = 0, -^ = 0, (i = l,2,...,n). (4.7) 
dt dx i0 

Notice again that X contains the unknown function W(t ,x ) dependent on i an d x , 
hence the derivatives are total derivatives. As the generalization of the case for ODE's, 
we demand that the RG/E equation should be compatible with the condition that the 
coordinate of the point of tangency becomes the parameter of the family of the surfaces;i.e., 

to = t, x = x. (4.8) 

Then the RG/E equation is now reduced to equations for the unknown function W, which 
will be shown to be the amplitude equations such as time-dependent Ginzburg-Landau 
equation. Here we remark that although Eq.(4.7) is a vector equation, the equation to 
appear below will be reduced to a scalar one; see subsection 6.2. It can be shown, at 
least for equations treated so far and here, that the resultant envelope functions satisfy 
the original equations uniformly up to 0(e p ); see also Ref.||. 



5 Simple examples 



In this section, we treat a few of simple examples of systems of ODE's to show the how 



the RG method works. The examples are the Duffing |K| equation of non- autonomous 
nature, the Lotka-Volterra||TTJ and the Lorenz[|l^] equation. The first one may be treated 
as a scalar equation. Actually, the equation is easier to calculate when treated as a scalar 
one. We give such a treatment in Appendix B. We shall work out to derive the time 
dependence of the solution to the Lotka-Volterra equation explicitly. The last one is an 
example with three degrees of freedom, which shows a bifurcation[TI|. We shall give the 
center manifolds to this equation around the first bifurcation of the Lorenz model. A 
general treatment for equations with a bifurcation will be treated in section 6. 
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5.1 Forced Duffing equation 

The forced Duffing equations are reduced to 



x + 2e7± + (f + ea)x + ehx 3 = ef cost, 
y + 2e^yy + (1 + ecr)y + ehy 3 = ef sin t. 



Defining a complex variable z = x + iy, one has 

eh 



z + 2e 7 i + (f + eo)z + — (3\z\ 2 z + z* 3 ) = efe u . 



We suppose that e is small. 

We convert the equation to the system 

d 



where 



( s -!„)»= -**•({,„;*) (J), 



u = ( „ ) > £ = 2 ' »7 = 



L 



V, 
f ' 

-f 0; ' 



and 



Let us first solve the equation in the perturbation theory by expanding 

u = u + eui + . . . , 

with Uj = t {ii,rji) (i — 0, 1, . . .). We only have to solve the following equat 
sively; 

(|-L„)u„ = 0, 

(|-io) Ul = -F(6,»it)(;), 

and so on. The solution of the zero-th order equation is found to be 

u (t;t ) = W(t )Ve l \ 
where U is an eigenvector belonging to an eigen value i of L , 

L U = fU, U=Q). 



The other eigenvector is given by the complex conjugate U*, which belongs to the other 
eigenvalue —i. We have made it explicit that the constant W may be dependent on the 
initial time to- m terms of the component, 

&(Mo) = W(t )e u , m (t;t ) = tW(t )e u . (5.10) 

Inserting these into F(£ , rj ; t), one has 

Ffo,w,t) = W(to)e u + ^W* 3 e- m , (5.11) 



with 



W(t ) = {a + 2i-f)W + —\W\ 2 W - f (5.12) 



We remark that the inhomogeneous term includes a term proportional to the zero-th order 
solution. Thus u x contains a resonance or a secular term as follows; 

Ui(Mo) = -^We i *{(t-t + ^)U-^U*}-^* 3 e- 3 ^(U-2U*). (5.13) 
In terms of the components 

6 (Mo) = ^W'(t-t ) + A W *V 3 *', 

?7i (Mo) = -^(t-to-^-^hW^e-™. (5.14) 

Adding the terms, we have 
u(t) ~ uo(t;*o) + eui(*;* ), 

= W(t )Ue u - e^We u {(t -t + ^)U - ^U*} - e^W* z e~ m {\J - 2U*), 
= u(t;t ). (5.15) 
In terms of the components, 

£(t;t ) ~ ^(t )e^ + e lw(t )e l '(t-to) + e^*V 3i ' = e, 

77(t;* ) ~ z^(t )e^-e^(t-t -2)-e^W*V 3i ' =77. (5.16) 

I lb 

Now let us construct the envelope ug(t) of the family of trajectories or curves u(t; to) — 
(£(t; to),fj(t; to)) which is parametrized with t ; u^(t) will be found to be an approximate 
solution to Eq. (5.3) in the global domain. According to section 2, the envelope may be 
obtained from the equation 

^ ;to) =0. (5.17) 



10 



In the usual procedure for constructing the envelopes, the above equation is used for 
obtaining t as a function of t, and the resulting t = t (t) is inserted in u(t; t ) to make 
the envelope function Ue(t) = u(t;t (t)). In our case, we are constructing the envelope 
around t = t , so we rather impose that 

t Q = t, (5.18) 

and Eq.(5.17) is used to obtain the initial value W(t ) as a function of to- That is, we 
have 



du(t; t ) 



dt 



o 



t =t 



dW TT it W itTT idW a l, TT T _ 3ehdW* „ 3it/TT 

= ^r Ue + e 27 e u + e 2^ e V u - u ) - i6^r e (u - 2U } - (5 - 19) 

Noting that the equation is consistent with dW/dt = 0(e), one has 
dW .e, A1/ v 

IT = ! 2 W(i) ' 

= .-{( ( x + 2. 7 )^(t) + T |^(t)| 2 ^(t)-/}. (5.20) 

This is the amplitude equation called Landau-Stuart equation, which may be also given 
by the RP method f| as a reduction of the dynamics. With this equation, the envelope 
trajectory is given by 

£ E (t) = W(t)e u + e^-W* 3 e- 3lt , 

lb 

VE(t) = i(W(t) + t\w[t))e lt -e%hW* 3 e~ 3a . (5.21) 

2 lb 



For completeness, let us examine the stationary solution of the Landau-Stuart equa- 
tion, briefly; 



W=(a + 2i-i)W + -eh\W\ 2 W - f = 0. 



(5.22) 



Writing W as 



W = Ae 



-i0 



(5.23) 



we have 



A 2 



(~hA 2 + a) 2 + 4 7 2 



(5.24) 



which describes the jumping phenomena of the Duffing oscillator. 
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5.2 Lotka-Volterra equation 



As another simple example, we take the Lotka-Volterra equation JIT 



x = ax — exy, y = —by + e'xy, (5.25) 

where the constants a, b, e and e' are assumed to be positive. It is well known that the 
equation has the conserved quantity, i.e., 

6 In \x\ + aln \y\ — (e'x + ey) = const.. (5.26) 



The fixed points are given by (x = 0,y = 0) and (x = b/e',y = a/e). Shifting and 
scaling the variables by 

x=(b + eO/e\ y = a/e + rj, (5.27) 
we get the reduced equation given by the system 

-L )u = -e^(_ 1 J, (5.28) 



d 
It 



L °=(°. n)- (5-29) 



where 

' // }' ~ \ a i) 

The eigen value equation 

L U = A U (5.30) 

has the solution 

A = ±iVab = ±iw, U = (^) • (5-31) 

Let us try to apply the perturbation theory to solve the equation by expanding the 
variable in a Taylor series of e; 

u = u + eui + e 2 u 2 H , (5.32) 

with Uj = t (^i,r]i). The lowest term satisfies the equation 

d 



(It -L )uo = 0, (5.33) 
which yields the solution 

u (t; t ) = W(t )e iuJt V + c.c, (5.34) 
12 



or 



fo = W(t )e^ + c.c, 77o = -^(iW^to)^ + c.c). 
Here we have supposed that the initial value W depends on the initial time t . 
Noting that 

( J^) — "U + c.c, 

with a — (1 — ib/uj)/2, one finds that the first order term satisfies the equation 



the solution to which is found to be 



iW 2 e 2iwt (aU + c.c.) + c.c. 



Ul= b 



iy 2 (aU+^-U*)e 2iw '+c.c. 



or 



,2b + iu„ r2 2iut 



3b 2 



(5.35) 



(5.36) 



(5.37) 



(5.38) 



(5.39) 



The second order equation now reads 

4 - L0V2 = t^t [{(6 - ico)\W\ 2 We iuJt + 3(6 + ^)W 3 e 3 ^} + c.c. 
at J 3b 2 L 



-1 / ' 



(5.40) 



We remark that the inhomogeneous term has a part proportional to the zero-th-order 
solution, which gives rise to a resonance. Hence the solution necessarily includes secular 
terms as follows; 



u 2 = 



b — iuj 



+ 



3b 2 

b + iuj 



W\ 2 w\a(t - t + i— )U + — trie*" 
I 2c; 2ia; J 



— H/ 3 (2aU + a*U*)e 3 ^' 



46 2 ?,CJ 

In terms of the components, one finds 
-i b 2 + uj 2 



+c.c. 



(5.41) 



6 = 

V2 = 



w 



6uj b 2 
\W\ 2 W 
66 3 
W 3 

+ 



W\ 2 W(t - t )e tut + -—{(3uj 2 - b 2 ) - 4i6cu}e : 



3iujt 



UJ^ 



+C.C. 



-(b 2 + uo 2 )(t - t ) + -{2buj + i(b 2 - uj 2 )} 



UJ 



Auit 



, o { 46 + — (36 w "I - c.c. 

06 Ll> 



(5.42) 



13 



The RG/E equation reads 



du 

wr ' (5 - 43) 



with to = t, which gives the equation for W(t) as 

~w = -^^ mw - (5 ' 44) 

If we define A(t) and 9(t) by W(t) = (A(t)/2i)expi9(t), the equation gives 

A(t) = const., 0(*) = -1— - (1 + — )w* + o , (5.45) 

24 

with being a constant. Owing to the prefactor % in r.h.s. of Eq. (5.44), the absolute 
value of the amplitude A becomes independent of t, while the phase 9 has a t-dependence. 
The envelope function is given by 



u * w = (Sw) = u((A 



(5.46) 

t =t,du/dt =0 



In terms of the components, one has 



£ E = AsmQ(t) - e^-(sin29(t) + ^cos26(t)) 

e 2 A 3 3u 2 - b 2 . , * 4o;6 

(sin 36 (t) - — = — -cos30(t)), 



32 u 2 b 2 v w 3u 2 -b 2 

c 2 A3 _ ,.,2\ ,2/13 



A -^^) cose(t) -T2^ sm0(t) 



T 

A 2 / 2b \ e 2 A 3 f 3b 2 — uj 2 X 



where 



2 y^2 ^2 

0(t)=(W + 0o, w = {l-i— + (5.48) 

24 or 

One sees that the angular frequency is shifted. 

We identify u E (t) = T]E(t)) as an approximate solution to Eq.(5.28). According 

to the basic theorem presented in section 4, Ue(i) is an approximate but uniformly valid 
solution to the equation up to 0(e 3 ). We remark that the resultant trajectory is closed 
in conformity with the conservation law given in Eq. (5.26). 

"Explicit solutions" of two-pieces of Lotka-Volterra equation were considered by Frame 
nj ; however, his main conceren was on extracting the period of the solutions in an average 



method. Comparing the Frame's method, the RG method is simpler, more transparent 
and explicit. The present author is not aware of any other work which gives an explicit 
form of the solution as given in Eq. (5.47,48). 
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5.3 The Lorenz model 



The Lorenz model [T^] for the thermal convection is given by 

£ = <r(-Z + r}), 

C = iv-K- (5.49) 

The steady states are give by 



(A) (£,77,0= (0,0,0), (B) (£, v , = (±yjb(r - 1), ±yjb(r - 1), r - 1). (5.50) 



The linear stability analysis Jl0| shows that the origin is stable for < r < 1 but 



unstable for r > 1, while the latter steady states (B) are stable for 1 < r < a(a + b + 
3)/(er — b — 1) = r c but unstable for r > r c . In this paper, we examine the non-linear 
stability around the origin for r ~ 1; we put 

r = 1 + n and /j = %e 2 , X — sg 11 /^- (5.51) 

We expand the quantities as Taylor series of e: 

u = ( 77 1 = eu! + e 2 u 2 + e 3 u 3 + • • • , (5.52) 

where Uj = *(^, rji, Q) (i = 1, 2, 3, . . .). The first order equation reads 

d 



dt 



where 







—0 


a 


1 


-1 









(5.53) 



(5.54) 



the eigenvalues of which are found to be 

Ai = 0, A 2 = -(j-l, A 3 = -b. (5.55) 
The respective eigenvectors are 

Ui=flj, U 2 =^-lj, U 3 =|^oj. (5.56) 

When we are interested in the asymptotic state as t — > 00, one may take the neutrally 
stable solution 

u 1 (t;t ) = W(t )XJ 1 , (5.57) 
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where we have made it explicit that the solution may depend on the initial time to, which 
is supposed to be close to t. In terms of the components, 



Zi(t) = w(t ), m(t)=w(t ), Ci(*) = o. 



(5.58) 



The second order equation now reads 



d_ 

dt 



- L ^j u 2 



t 

-£iG I = iv-u,. 



which yields 



W 2 

u 2 (t) = —U 3 , 



or in terms of the components 



W 2 

6 = V2 = 0, C2 = — • 



Then the third order equation is given by 

( ° ^ 

-x6 - 6C1 - 6C2 



(5.59) 



(5.60) 



(5.61) 



-±—{ x W- ^(aUi-Ua), 

1 + <7 



(5.62) 



which yields 



u 3 = t^-(xW - V)M* - to + t^-)Ui - T"^U 2 } 
1 + a b l + o" 1 + cr 



(5.63) 



Thus gathering all the terms, one has 



u(t;fo) = eW(to)Ui + jW{t ) 2 U 3 



+ -rr^(xW(t ) - l W(t Q f ){a(t - t + - -^U 2 },(5.64) 

1 + cr 1 + a l + o 



up to 0(e 4 ). The RG/E equation now reads 



= 



du 

dt 



dW, 



t =t 



.t 2 „dW, 



= t ^—\J l + 2-W^—V,--^—e\ X W- V 3 )Ux, 
dt b dt 1 + cr b 



(5.65) 



up to 0(e 4 ). Noting that one may self-consistently assume that dW/dt = 0(e 2 ), we have 
the amplitude equation 



(5.66) 
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With this W(t), the envelope function is given by 



u E (t) = u(t;t = t), 



eW(t)Vi + ^W(t) 2 V 3 + f , 2 (xW(t) - ]-W(tf)(aV 1 - U 2 ),(5.67) 
o 1 + ar o 



or 



^(t) = eW(t), 

VE (t) = eW{t) + -^-{ X W{t)-\w{tf), 
1 + a o 

( E (t) = jW(tf. (5.68) 

We may identify the envelope functions thus constructed as a global solution to the Lorenz 
model; according to the general theorem given in section 4, the envelope functions satisfy 
Eq.(5.49) approximately but uniformly for Vt up to 0(e 4 ). 

A remark is in order here; Eq.(5.68) shows that the slow manifold which may be 
identified with a center manifold |TD[ is given by 



1 + a b(l+ a) b 

Notice here that the RG method is also a powefull tool to extract center manifolds in 
a concrete form. It is worth mentioning that since the RG method utilizes neutrally 
stable solutions as the unperturbed ones, it is rather natural that the RG method can 
extract center manifolds when exist. The applicability of the RG method was discussed 
in 0| using a generic model having a center manifold, although the relation between the 
exitence of center manifolds and neutrally stable solutions is not so transparent in their 
general approach. 



6 Bifurcation Theory 



In this section, we take generic equations with a bifurcation. We shall derive the Landau- 
Stuart and Ginzburg-Landau equations in the RG method. In this section, we shall 
follow Kuramoto's monograph || for notations to clarify the correspondence between the 
RG method and the reductive perturbation (RP) method. 



6.1 Landau-Stuart equation 



We start with the n-dimensional equation 



— = F(X,t;/i). (6.1) 
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Let X (/i) is a steady solution 

F(X (/i);/i) = 0. 
Shifting the variable as X = X + u, we have a Taylor series 

du 

— = Lu + Muu + A^uuu + • • • , 

at 

where we have used the diadic and triadic not at ions M: 



dX 3 



1 d 2 F 

x=x ,- fc ZoXjdXk 



(iVuuu), = X! 



1 <9 3 F,- 



^dXjdX^Xt 



X=Xn 



UjU k Ui. 



(6.2) 



(6.3) 



X=Xn 



(6.4) 



We suppose that when \i < 0, X is stable for sufficiently small perturbations, while 
when n > 0, otherwise. We also confine ourselves to the case where a Hopf bifurcation 
occurs. We expand L, M and N as 



L = L + fiL 1 + ---, M = M + fiM 1 + ---, N = N + fiN 1 + 
The eigenvalues X a (a = 1, 2, . . . , n) of L are also expanded as 

A° = A° + //A? + ■ • • , 

with 

A)U a = AgU Q . 



(6.5) 



(6.6) 



(6.7) 



We assume that Aj = — Ag are pure imaginary, i.e., Aj = iu , and ReAg < for a 



3,4,.... 



Defining e and x by e = J and x = sgn/x, we expand as 



u = eui + e 2 u 2 + e 3 u 3 + 



(6.8) 



The Uj (i = 1, 2, 3, ...) satisfies 



d_ 

dt 
d 



(di~ L °) U2 = M ° UlUl ' 

-L )u 3 = xLiUi + 2M Uiu 2 + iV uiUiUi, 



dt 



(6.9) 



etc. 
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To see the asymptotic behavior as t — > oo, we take the neutrally stable solution as the 
lowest one around t = t ; 

ui(t; t ) = W(t )Ue wot + c.c., (6.10) 

where c.c. stands for the complex conjugate. With this choice, we have only two degrees 
of freedom for the initial value W(t ). 

The second order equation is solved easily to yield 

u 2 (Mo) = (V + ^(t ) 2 e 2 ^ + c.c.)+V |^(to)| 2 , (6.11) 

where 

V + = -(L - 2iu Q )- 1 M UU ', V = -2Lo 1 M UtJ, (6.12) 

with U being the complex conjugate of U.Q Inserting u x and u 2 into the r.h.s of Eq. 
(6.9), we get 

L )u 3 = {x^WU + (2M \JV + + 3N U\JXJ)\W\ 2 W}e iuJot + c.c. +h.h., 



dt 



Ae^ ' + c.c. + h.h., (6.13) 



where h.h. stands for higher harmonics. So far, the discussion is a simple perturbation 
theory and has proceeded in the same way as given in the RP method except for not 
having introduced multiple times. 

Now we expand A by the eigenvectors U Q of L as 

A = 5> a U a , (6.14) 



where 



Here U* satisfies 



A a = U^A. (6.15) 



U;L = XqLq, (6.16) 

and is normalized as U^Uq, = 1. 



Then we get for 113 

u 3 (t; t ) = {Ax(t -t + 5)\J+J2 . A \ a V a }e luJot + c.c. + h.h.. (6.17) 

a _^ 1 iloq - A 



7 In other sections, we use the notation a* for the complex conjugate of a. In this section, * is used 
for a different meaning, following ref. [0; see Eq. (6.16). 



19 



The constant 5 is chosen so that the coefficient of the secular term of the first component 
vanishes at t — to. Note the appearance of the secular term which was to be avoided in 
the RP method: The condition for the secular terms to vanish is called the solvability 
condition which plays the central role in the RP method[|J. 

Thus we finally get 

u(t; t ) = {eW(t Q )U + e 3 ^ - t + S)U + £ , Aa Q U a )}e^ f + c.c. + h.h.. (6.18) 



The RG/E equation 



du 

dt 



= 0, (6.19) 

t =t 



yields 

dW 



e 2 A 1 , 



dt 

= e 2 [xU*LxUiy + {2U*M Q UV + + 3U*N UU\J}\W\ 2 W], (6.20) 

up to 0(e 3 ). Here note that the terms coming from h.h. do not contribute to this order 
because dW/dt is 0(e 2 ). The resultant equation is so called the Landau-Stuart equation 
and coincides with the result derived in the RP methodlBI. 



6.2 The Ginzburg-Landau equation 

We add the diffusion term to Eq.(6.1); 

^ = F(X) + DV 2 X, (6.21) 
where D is a diagonal matrix. Let Xo be a uniform and steady solution. 
Shifting the variable X = X + u as before, we have 

y = Lu + Muu + iVuuu H , (6.22) 

with 

L = L + DV 2 . (6.23) 



Then using the same expansion as before, we have the same equation for ui,U2 and 
u 3 as given in Eq.(6.9) with L being replaced with L = L + W 2 . 
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To see the asymptotic behavior as t — > oo, we take the neutrally stable uniform solution 
as the lowest one around t = t and r = r ; 

ui(*, r; t , t ) = W(t , r )Ue^* + c.c. (6.24) 

With this choice, we have only two degrees of freedom for the initial value W(t , r ). 

The second order equation is solved easily to yield the same form as that given in 
Eq.(6.11). Inserting u x and u 2 into the r.h.s of Eq. (6.9) with L Q replaced with L , we 
have 





Of 



-L Ju 3 = {x^iW / U + (2M UV + + 3iVoUUU)|W / | 2 iy}e^ ot + c.c. + h.h., 

= Ae iuJot + c.c. + h.h.. (6.25) 
Then we get for 113 in the spatially 1-dimensional case, 



u 3 (t;*o) 



A 1 {a 1 (t -to + 6)- f zr V -4 + 0>U + E — ^ 



+c.c. + h.h., 



(6.26) 



with c\ + C2 = 1. We have introduced constants 5 and 5' so that the secular terms of the 
first component of U3 vanish at t — to and x = xq. Note the appearance of the secular 
terms both t- and x-directions; these terms were to be avoided in the RP method with 
the use of the solvability condition. 



Adding all the terms, we finally get 
u(Mo) = 



eW(t , x )U + e 3 {A 1 ( Cl (t - t + 5) - -jD~\x 2 - x\ + <f'))u 

+ E -A^ u «} 

a ^ x i^o - A 



e iuJot + c.c. + h.h. 



(6.27) 



up to 0(e 4 ). The RG/E equation^ 



du 

dt 



t =t 



du 

dxo 



0. 



(6.28) 



Xq = X 



yields 



e 2 c l A l + 0(e 3 ), D- 



dt 



dx 



-e 2 arc 2 A 1 + 0(e 3 ). 



(6.29) 



We remark that the seemingly vector equation is reduced to a scalar one. Differentiating 
the second equation once again, we have 



D 



d 2 W 
dx 2 



-e 2 c 2 A 1 + 0{e 3 ). 



(6.30) 



^See section 3. 
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Here we have utilized the fact that dW/dx = 0(e 2 ). Noting that c\ + C2 = 1, we finally 
reach 



m_ D cPw_ 2A 

dt dx 2 1 ' 

= e 2 [xU*L x \JW + {2U*M UV+ + 3U*N UXJ\J}\W\ 2 W}, (6.31) 

up to 0(e 3 ). This is so called the time-dependent Ginzburg-Landau (TDGL) equation 
and coincides with the amplitude equation derived in the RP method||. 

We have seen that the RG method can reduce the dynamics of a class of non-linear 
equations as the RP method can. Therefore it is needless to say that our method can be 
applied to the Brusselators|15j, for instance, and leads to the same amplitude equations 
as the RP method|| does[@. 



7 A brief summary and concluding remarks 

In this paper, we have shown that the RG method of Goldenfeld, Oono and their collab- 
orators can be equally applied to vector equations, i.e., systems of ODE's and PDE's, as 
to scalar equations. || |], |], || We have formulated the method on the basis of the classical 
thoery of envelopes, thereby completed the argument given in [^, |6|] . We have worked out 
for some examples of systems of ODE's, i.e., the forced Duffing, the Lotka-Volterra and 
the Lorenz equation. It has been also shown in a generic way that the method applied to 
equations with a bifurcation leads to the amplitude equations, such as the Landau-Stuart 
and the (time-dependent) Ginzburg-Landau equation. 

Then how about the phase equations §?: The phase equations describe another re- 
duced dynamics. The basis of the reduction of the dynamics by the phase equations lies in 
the fact that when a symmetry is broken, there appears a slow motion which is a classical 
counter part of the Nambu-Goldstone boson in quantum field theory. We believe that 
if the phase equations are related to slow motions of the system at all, the RG method 
should also leads to the phase equations. It is an interesting task to be done to show that 
it is the case. 

There is another class of dynamics than those described by differential equations, i.e., 
difference equations or discrete maps. It is interesting that a natural extension of the 
RG/E equation to difference equations leads to a reduction of the dynamics. |9|] This fact 
suggests that the RG method pioneered by Goldenfeld , Oono and their collaborators 
provides one of the most promising candidate for a general theory of the reduction of 
dynamics, although it is certain that such a mechanical and general method is often tedious 
in the actual calculations.^ As an application of the reduction of difference equations, it 

9 It should be mentioned that there are other methods jLTj , |l8| for the dynamical reduction as promising 
as the RG and RP method are. 
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will be interesting to see whether the coupled map lattice equations as systems of non- 
linear difference equations fl9| can be reduced to simpler equations by the RG method. 
We hope that we can report about it in the near future. 
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Appendix A 



In this Appendix, we give a quick review of Goldenfeld, Oono and their collaborators' 
prescription for the RG method. Then we summarize the problems to which a mthematical 
reasoning is needed in the author's point of view. 

We take the following simplest example to show their prescription: 

d x dx , . 

- 2 + e- + x = 0, (A.l) 

where e is supposed to be small. The exact solution reads 



x(t) = Aexp(~t) sin(y l-jt + 6), (A.2) 
where A and 9 are constant to be determined by an initial condition. 

Now, let us blindly apply the perturbation theory expanding x as 

x{t) = x (t) + ezi(t) + e 2 x 2 {t) + .... (A.3) 
The result is found to be|| 

x(t;t ) = A sm(t + 6 )-e^(t-t )sm(t + 6 ) 

+e 2 ^{(t - t ) 2 sin(t + 6 ) - (t - t Q ) cos(t + 9 )} + 0(e 3 ). (A.4) 

o 

Now here come the crucial steps of the Goldenfeld et al's prescription: 

(i) First they introduce a dummy time r which is close to t, and "renormalize" x(t;to) 

by writing t — t = t — t + t — t ; 

x{t,r) = A{r)sm{t + 6{T)) - e ^p-(t - r)sm(t + 6(t)) 

+e 2^1{( t _ T f sin ( t + o( T )) -(t-r) cos(t + 9(r))} + 0(e 3 fA.5) 

o 

with 

x(r,r) = i4(r)sin(r + 0(r)). (A.6) 
Here Aq and 8q have been multiplicatively renormalized to A(t) and 0(t). 

(ii) They observe that r is an arbitrary constant introduced by hand, thus they claim 
that the solution x(t, r) should not depend on r; namely, x(t, r) should satisfy the 
equation 

= 0. (A.7) 

dr 

This is similar to the RG equation in the field theory where r corresponds to the 
renormalization point r; hence the name of the RG method. 
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(iii) Finally they impose another important but a mysterious condition that 



T = t. 



(A.8) 



From (ii) and (iii), one has 




(A.9) 



which gives 



A(t) = Ae~ €T/2 , 



9(t) 




(A.10) 



where A and 9 are constant numbers. Thus, rewriting r to t in x(t), one gets 



x(t,t) = Aexp(--t) sin((l - —)t + 6). 



(A.11) 



They identify x(t,t) with the desired solution Then one finds that the resultant 

x(t) is an approximate but uniformly valid solution to Eq.(A.l). In short, the solution 
obtained in the perturbation theory with the local nature has been "improved" by the 
RG equation Eq.(A.7) to become a global solution. 

But what have we done mathematically? what is a mathematical meaning of the 
" renormalization" replacing t with the extra dummy time r? Can't we avoid the "renor- 
malization" procedure to solve a purely mathematical problem? Why can we identify 
x(t,t) with the desired solution?; with r being a constant, x(t,r) can be a(n) (approxi- 
mate) solution to Eq. (A.l), can't it? In other words, when the operator d/dt hits the 
second argument of x(t,t), what happens? 

In Ref. 0, it was shown that the "renormalization" procedure to introduce the extra 
dummy time r is not necessary. Furthermore, it was clarified that the conditions (ii) and 
(iii) are the ones to construct the envelope of the family of the local solutions obtained in 
the perturbation theory; x(t; t) is the envelope function of the family of curves given by 
x(t;to) where to parametrizes the curves in the family. Furthermore, it was shown that 
the envelope function x(t, t) satisfies the orginal equations approximately but uniformly; 
the hitting of d/dt on the second argument of x(t,t) does not harm anything. In short, 
the prescription given by Goldenfeld, Oono and their collaborators is not incorrect, but 
the reasoning for the prescription is given in |^, |5| and will be more refined in the present 
paper. In Ref. 0, a simplification of the prescription and its mathematical foundation is 
given for PDE's. 
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Appendix B 

In this Appendix, we solve the forced Duffing equation without converting it to a 
system. It is easier to solve it in this way than in the way shown in the text. 

We start with Eq. (2.6) 

z + 2e 7 i + (1 + ea)z + y (3\z\ 2 z + z* 3 ) = efe l \ (B.l) 

where e is small. 
Expanding z as 

z = z + ezx + t 2 z 2 + ■■■, (B.2) 
one gets for z in the perturbation theory 

Q /i 1 

z(t; t ) = W{t )e* + e(t - t ){f - W{a + 2z 7 ) - —\W\ 2 W}e u + e -W*(t) 3 e m + 0(e 2 Y.B.3) 

2 16 

Note that there exists a secular term in the first order term. 
The RG/E equation reads|5|] 

with to — which leads to 

W = -e{a + 2i 1 )W- ^eh\W\ 2 W + ef (B.5) 

up to 0(e 2 ). Here we have discarded terms such as edW/dt, which is 0(e 2 ) because 
dW/dt = 0(e). The resultant equation for the amplitude is the Landau-Stuart equation 
for the Duffing equation. The envelope is given 

z E (t) = z(t; t = t) = W(t)e lt + ^-W* 3 e 3U + 0(e 2 ). (B.6) 

16 

We identify Zs{t) with a global solution of Eq.(B.2), and x(t) = Re^] and y(t) = Im[^ s ] 
are solutions to Eq.(B.l). As shown in the text, Vt, Zs(t) satisfies Eq.(B.2) uniformly up 
toOfe 2 ). 
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